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1.  INTRODUCTION 


-^?The  use  of  lidar  to  infer  information  about  battlefield  obscurants 
such  as  smoke  and  dust  has  received  considerable  attention  in  recent  years. 


Review  articles,^-’ ^/feasibility  studies and  experiments'*’^*  have  demon- 

_ ' 

strated  thaf^state-of-the-art  lidars  can  provide  relative  transmission  and 

* 

spatial  concentration  data  for  three  dimensional  obscurant  clouds.  In 
spite  of  this,  the  theoretical  models  necessary  for  interpreting  the  lidar 
return  have  not  contained  dependencies  on  all  of  the  relevant  transmitter 
and  receiver  parameters,  and  have  been  almost  entirely  limited  to  fi^st 
order  scatter.  Since  the  conditions  for  these  models  to  be  valid  are  not 
always  met,  there  appeared  to  be  a  need  for  a  more  fundamental  theoretical 
basis  which  would  reproduce  the  usual  equations  as  an  approximation  and 
would  give  a  clear  basis  for  understanding  their  limitations  and  general¬ 
izing  them  as  required.  Providing  this  model  was  one  of  the  primary  aims 
of  the  contract. 

In  addition,  analyses  of  the  inversion  problem  have  been  largely 
qualitative  and  lacking  in  explicit  mathematical  algorithms  for  deriving 
information  about  the  scattering  medium  from  the  lidar  return.  The  work 
of  Lamberts^  on  finding  statistical  correlations  between  the  lidar  return 
and  independently  observed  aerosol  scattering  data  is  an  exception,  but  it 
is  unclear  how  a  quasi-real  time  algorithm  could  be  developed  along  these 
lines.  For  this  reason  another  important  aim  of  the  contract  was  to  de¬ 
rive  explicit  algorithms  using  the  lidar  signal  either  by  itself  or  in 
conjunction  with  supporting  measurements  to  derive  information  about  the 
scattering  medium. 

It  is  essential  to  specify  the  type  cf  information  required  about  the 
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obscuring  medium;  for  tactical  applications  usually  only  the  optical  trans¬ 
mission  is  important.  In  this  case  it  appears  possible  to  perform  a  re¬ 
latively  simple  integration  of  the  lidar  return.  For  field  studies  on  the 
kinetics  of  explosive  rounds,  however,  a  detailed  space-time  history  of 
the  debris  concentration  is  needed  and  perhaps  even  an  estimate  of  particle 
size  distributions  and  relative  numbers  of  scattering  constituents.  These 
diagnostic  goals  are  much  harder  to  implement  and  either  supporting  local 
measurements  of  the  particle  scattering  cross  sections  or  use  of  multiple 
wavelengths  is  required  in  addition  to  the  lidar  signal.  Part  of  the 
contract  involved  attempting  to  specify  clearly  what  additional  informa¬ 
tion  would  be  needed  to  perform  these  additional  tasks. 

In  section  2.1  a  general  solution  to  the  radiation  transport  equation 
is  obtained  as  an  expansion  in  the  ratio  of  backseat ter ing  to  extinction 
coefficients.  This  is  valid  because  computer  calculations  of  Mie  scattering 
for  particles  with  refractive  indices  such  as  smoke  and  dust  show  this 
ratio  is  much  less  than  one;  thus  the  series  converges  quickly. 

In  section  2.2  approximations  of  small  angle  scattering  and  Gaussian 
aperture  weighting  functions  are  combined  to  yield  the  dependence  of  the 
first  order  backscattered  power  on  the  lidar  receiver  aperture  size  and 
field-of-view,  as  a  function  of  time. 

In  section  2.3,  assuming  a  Gaussian  transmitter  beam  and  a  three  di¬ 
mensional  density  distribution  for  the  smoke  cloud,  the  expression  developed 
in  section  2.2  is  integrated  to  yield  an  analytical  expression  for  in 
terms  of  the  parameters  of  the  lidar  system  and  the  environmental  parameters 
of  the  smoke  cloud.  Wt  then  introduce  a  model  for  the  return  from  a  target 
positioned  beyond  the  smoke  cloud  and  obtain  the  total  lidar  signal  for 
the  cloud-target  combination.  This  latter  model  is  important  for  assessing 
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the  ability  of  the  smoke  cloud  to  obscure  a  target. 

Following  an  outline  of  our  inversion  approach  in  section  3.1  we  de¬ 
rive  an  analytic  expression  for  the  path  transmission  as  a  function  of 
the  lidar  signal,  system  constants,  and  backscatter-to-extincticn  ratio 
a  in  section  3.2. 

The  sensitivity  of  the  transmission  expression  to  noise  and  uncer¬ 
tainties  in  a  is  investigated  in  section  3.3  using  a  combination  of  analy- 
tical  and  simulation  techniques.  In  section  3.4  we  consider  means  of  en¬ 
hancing  the  inversion  algorithm  performance  by  additional  measurements  in 
a  small  region  of  the  propagation  path  of  scattering  parameters.  Finally, 

t 

we  outline  an  analytical  approach  for  using  multiple  wavelengths  to  infer 
particle  size  distributions  and  concentrations  in  section  3.5. 

Section  4  summarizes  our  conclusions  of  how  well  the  proposed  inver- 
sion  techniques  would  perform  and  suggests  additional  work  to  improve  the 
multiple  wavelength  approach. 
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2>  MULTIPLE  SCATTERING  THROUGH  A  SPATIALLY  VARYING  TURBID  MEDIUM 

2.1  General  Formalism 

In  this  section  we  present  a  formal  solution  to  the  radiative  trans¬ 
port  equation  tor  a  spa'  tally  varying  medium.  Although  the  application 
to  lidar  involves  predominantly  single  scattered  radiation  for  smoke  and 
dust  aerosols,  it  is  useful  to  base  the  calculations  on  a  more  general 
formulation  that  is  capable  of  straightforward  extention  to  other  situ¬ 
ations. 

The  radiative  transport  equation  for  a  general  scattering  medium  has 
the  form 


1  ?J 
c 


(r,n,t)  +  n  •  J(r,n,t)  =  -e(r,t)  J(r,n,t) 

0(r,t)  r 


4ir 


J  dun,  x(n*n 

4ir 


’)  J(£»n' ,t)  f  Sq  (r,n,t)  ,  (1) 


where  J  is  the  power  radiated  into  direction  n_  at  location  £  at  time  t, 
e  is  the  extinction  coefficient,  0  is  the  scattering  coefficient  and  Sq 
is  a  source  term.  The  solid  angle  integration  is  carried  out  over  the 

scattering  phase  function  x  which  represents  the  relative  power  scattered 

,  8 
into  direction  n*ii  .  To  solve  Eq.  (1)  define 


a(r,t) 


3(r,t) 


e(r,t)  - 


<  1 


and  set 


J0r,n,t)  =  ^  [a(r,t)]k  Jk<£,n,t} 


(2) 


k-0 


I 
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Substitution  into  Eq.  (x)  produces  the  coupled  set  of  equations 


l3Jk 


37"  +  Jk  +  £  (£*.£.)  Jk  =  Sk(r,n,t) 


where 


sk^l’-,c^  “  »  k  =  0 


e(r,t) 


f  X  (£•£.')  Jk_1(l.»n'»t)  »  k 


>  0  . 


These  equations  can  be  solved  iteratively  using  the  input  source  function 


Sn.  An  equivalent  procedure  is  to  look  for  solutions  of  the  form 


Jk(r,n,t)  -  c  j  At'  J  d3  J  d“ni  G0(r,n,t;  r',n’,t')  Sk(r',n'  ,t')  , 


where  G0  is  the  Green's  function  or  propagator  for  scattered  radiation. 


Using  the  method  of  characteristics  Gq  can  be  shown  to  be 


c 

GQ(r,n,t;  _r ' ,  _n '  ,  t')  =  exp£-  c.  J* Ax  e  (r-nc ( t-r ) )j  63(;r-£'  -nc(t-t'))  6^(ri-n' ) 


where  we  have  assumed  e  to  be  time  independent  (or  slowly  fluctuating  over 


the  radiation  transit  time).  We  now  define  Gk>  the  Green's  functi  >n  for 


k  order  scattered  radiation,  as 


■  c 


f At'  /  d3  r'  j  dcun,  &k(r,n,t;  r'.n'.t')  S0(r',n',t') 

-»  ^  An  ,, 


An  expression  for  Gk  can  be  found  by  substituting  the  definition  of  Sk 


into  Eq.  (4)  and  applying  Eq.  (6);  the  result  is 


msssKssmsm 


r/idfe 


V>J'-“-.--»10'.(  1U  ■.•*»* r-r,  ‘.**v 


c  •  .E.'  >£.'  ,  t')  =  c 


idt"F r". 


d“n..  GQ(r,n,t;  r",n",t'') 


s(r") 


4tt 


f dcon„t  x  (n"*n"'  )  G^  (r",nm  ,t";  r',n'}t')  . 


•/ 

Att 


(7) 


Thus,  can  be  imputed  as  an  integral  over  the  arbitrary  source  distri¬ 

bution  of  as  given  in  Eq.  (7)  iu  terms  of  lower  order  G^'s.  Since 
the  G's  are  independent  of  the  source  dis^j-L  *tion,  the  problem  has  for¬ 
mally  been  solved  for  arbitrary  S  q. 

As  an  example  of  using  Eq.  (7)  we  find  for 

t 

G^r.n.t;  r',n',t')  -  “X<n*n')  f  dt"  T^t.t")  e  (r-nc(t-t")) 

' 

x  ^(t'^t^Ar-r'-ncN")  -n'c(t"-t')),  f8) 


where 


Tj,  (t , t")  =  exp 


t 

[-c  J"  dx  e  (r-nc(t-x))j 

c11 

T2(t",t'  )=  exp£-  c  J  dx  e  (£.'+2.'c(x-t' ))j  . 

Equation  (8)  can  be  simplified  substantially  in  the  case  of  a  uniform 
medium*,  namely 


uniform 


G1(r,n,t;  j:'  ,n'  ,t' ) 


medium 


c  e 


X  (£•£’)  _cc(t_t.) 


An: 


J  dt"  5  3(r-r’  -  n< 

V 


-r'  -  nc(t-t")  -n'cCf'-t’)) 


X(n.n')  *  <**’  ^1  c(t-Co>  i'c(t0“t* )) 

ce  - — —  e  -  — 


4ir 


ln3  -n^l 


(9) 
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where 


z-z'  -  n3  ct  +  n3'  ct' 
c(n3’  -  n3) 


and 


L  =  (£.»  '•  /  , 


n 


(*_!_,  n3) 
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•  2.2  Application  to  Backscatter 

Equations  6  and  7  allow  the  radiance  to  be  computed  formally  by 
orders  of  scattering  for  arbitrary  source  functions  and  spatial  distri- 

^  but  ions  of  scatterers.  In  practice  simplifying  assumptions  are  needed 

to  reduce  the  complexity  of  the  required  integrations.  In  this  section 
we  construct  an  expression  for  the  first  order  backscatter  power  P^(t) 

#  by  assuming  1)  the  small  angle  approximation  and  2)  the  collection 
aperture  and  bistatic  offset  in  the  receiving  plane  as  well  as  the  field- 
of-view  can  be  approximated  as  Gaussian  functions. 

•  Specializing  Eq.  8  to  the  backscatter  case  (z  **  z*  =*  0)  and  per¬ 
forming  the  t"  integration  gives 

c  X(n‘n’) 

*  G1  ( P>Xl* 1 5  p_' , n ’ , t * )  - e(r  -  nc(t-tQ))  T^tjtg)  T^t^t') 

l»3  “  n3l 

x  (p_  -  p_'  -  n^c (t  -  tQ)  ~  nj_c(t0  -  t’))  (10) 


where 


n3tf 


n 


T 

3 


"‘S  • 


For  the  applications  in  mind  the  lidar  transmitter  beam  divergence  and 
receiver  field-of-view  angles  are  small  ('v  milliradians) .  This  justifies 
making  the  small  angle  approximation  n  j_  -  £  >  -  £' .  This  gives 

tQ  «  (t  +  t’)/2. 

Also, 


m 


V 


2  -<L-b)2/4 

P  e 


Substitution  of  Eq.  11  gives 


f.\ 


-02/e2 

e  R  G1  (£.>£.» P.' ,n' ,t’).  (14) 


Pl(a',e',f)  -  f  ^  s^'  +  9*  c<4-^-)[T2(to,t')]2 

«  /j2ie  '  R  exP  [-(e-o  +  1  c<-"i"t’)J/4] 


where 


E0  5  £■'  +  i 


+  .'  c<i-=^--b  . 


Performing  the  integral  gives 


>>!<£.' .±V>  -  !  ^  6(e’  +  I'  [V'o-^]2 


*4  4 

4  *  4  (0<H^y 


exp - 


4  -  4<h^)2 


.(16) 


Equation  13  with  16  produces  the  desired  expression  for  ?1(t)j 

t 

piU)  -  f-  x(-D  a2  e2  fit’  /dV  /dV  s0(p,»l,*t,)  B(e.’  +  I*  c<* t^"5) 


[T2(t0,t’)l2  r  le.’  +0'  c<^-)  -b|2 

- —  eXp - - - 5 — 

(t  -  t ' )1  2  ,  2  r  (t  -  t'Vr 

C  2  j  aR  +  0R  LC  2  J 
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2 . 3  Coherent  transmitter  example 

We  apply  the  first  order  power  expression  Eq.  17  to  the  case  of  a 
narrow  coherent  Gaussian  transmitter  beam.  As  examples  we  construct 
analytic  expressions  for  the  power  collected  from  a  uniform  medium,  and 
present  numerical  results  for  a  three  dimensional  Gaussian  cloud. 

Assume  for  the  source  function 


^e'^Vd'  +EVf) 

Tra- 


ClS) 


where  E(t’)  specifies  the  temporal  shape  of  the  pulse.  Equation  18  re¬ 
presents  a  coherent  Gaussian  transmitter  wavefront  with  1/e  amplitude  a^ 
and  focal  length  f.  Substitution  of  Eq.  18  into  17  gives 


2  2 

c2  aR  6R 

VC)  "  8  2 


rit. - mn -  r 

l  J 


,2.  2 

.2  ,  "p  ^aT 

d  e 


g^o'(l  -  c(t2~  fc,))  +  e3  c(t  ~  ■t,')^T(t,t,)j 


exp 


k.(l  .  ^p)  .  b! 


where 


Tvt.t')  ■  exp 


2  . 

aR  +  9 


t+t! 

2 


2  Let  -  f)l2 
R  [  2  J 


(19) 


-c  j"  dt  e^p_’^l  "  c(x  ~  -^)-f  e3  c (t  -  t')^ 


L  f 


(20) 


( 
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In  general,  for  arbitrary  e(r),  Eq.  19  must  be  evaluated  numerically. 
For  the  special  case  e(r)  =  e  ,  a  constant,  P^(t)  becomes 


p^t)  -  x(-i)a^e 


R  3  fd 

R  0  J 


,  E(t')  e 


-e  c(t  -  t') 


2  2  1 

SR  +  0r[ 


[‘WJ' 


,  2  2 
+  aTY 


4 ♦  4  -  4v2 


where 


_  ,  c(t  -  t’)  ,  „  _ 

Y  =  1 - -  and  =  a  eq 


By  the  changes  of  variables  £'  =  c(t  -  t’),  £  =ct,  Eq.  21  may  be  written 


as  a  convolution 


?x(0  -  JdV 


E(£  -  £’)  P16(C) 


where 


P1S<« 


X(-l) 


a2  e2  6  e'E°5 

R  R  Q  _ 


4 +  4  (f>  +  4 (1  - 1?)2 


x  exp  |  - 


4 +  4. 4>  +4(1 


is  the  impulse  response  corresponding  to  a  delta  pulse  in  tine.  The  two 
limiting  cases  of  5  +  0  and  £  >>  f,  a^  are 


i 


13 


P1S<°> 


C  *(-l>46R 

8  2  2 

a_  +  a_ 

R  T 


6n  exp 


2  ,  2 
a„  +  am 


(24) 


11m 
£  >>  f,a. 


P1S<5) 


R 


v 

2 

2  2  £ 

K+6?<f> 


-eQC 


exp 


(25) 


( 


( 


r 


c 


c 


€ 


C 


with  0_  "  -a„/ f.  Notice  that  Eq.  25  reduces  to  the  usual  lidar  equation 

£nc  l 

as  b  +  0.  Figure  1  plots  log  [P^(£)  e  /-g-  cx(-1)Bq]  versus  £  =  ct  for 

_3 

various  values  of  the  offset  b  using  aD  =  10  cm,  a_  -  1  cm,  0O  -  0_  *  10  . 

K  T  _  K  T 

Vf  2  2 

0  —  +  0  ;  for  longer 
it  T 

2 

£  the  1/ £  term  dominates. 

Although  Equations  22  and  23  provide  insight  into  the  dependence  of 
lidar  signal  on  system  parameters,  they  are  not  adequate  to  model  the 
return  from  localized  distributions  of  scatterers.  For  this  reason  a 
computer  program  was  developed  to  perform  the  integrations  in  Eq.  19  for 
arbitrary  pulse  shapes  and  smoke/dust  distributions.  In  addition,  it  is 
often  useful  to  estimate  the  ability  of  a  lidar  to  detect  a  target  through 
an  obscuring  cloud.  For  simplicity  we  model  an  infinite  Lambertian  plate 
target  with  arbitrary  tilt  with  respect  to  the  transmitter  beam  axis. 
Appendix  3  derives  the  power  returned  from  this  target,  and  Figure  2 
plots  the  relative  received  power  versus  time  for  a  20  nsec  square  pulse 
incident  on  the  target  at  various  tilt  angles. 

Figure  3  plots  the  returns  from  a  10.6  pm  lidar  of  the  target  at 
normal  incidence  at  1  km  range  illuminated  through  a  three  dimensional 
Gaussian  smo’/c  cloud  centered  at  Zq  ■  .5  km  and  having  1/e  radius  a  ■  50  m. 


Figure  2.  Received  power  from  tilted  Lambertian  target 
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The  extinction  coefficient  e(r)  is  taken  to  be 


2,  2. 


e(r)  =  ea  +  eQ  exp  [-  -  r Q)/a  ] 


(26) 


where  »  [0,  0,  z^]  and  f.&  is  the  ambient  extinction.  Table  1  summarizes 
the  system  and  scattering  medium  parameters  U3ed  in  the  computer  calculation. 
The  different  plots  in  Figure  3  represent  different  values  of  Gq,  the  peak 
extinction  at  the  cloud  center.  Increasing  Eq  has  the  effect  of  reducing 
the  target  return  and  shifting  the  peak  of  the  cloud  return  toward  near 
ranges.  This  is  due  to  the  increased  attenuation  of  the  pulse  as  it  pene¬ 
trates  the  cloud  center.  For  the  case  of  largest  Eq  shown  here  the  apparent 
cloud  center  is  shifted  from  500  m  to  about  435  m.  This  effect  could  have 
an  important  bearing  on  using  lidar  to  map  cloud  concentrations. 
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Table  1 

PARAMETERS  USED  IN  LIDAS  SAMPLE  CALCULATION 

Lidar  System 


Symbol 

Meaning 

Value 

P^t) 

1st  order  Lidar  power 

Computed 

• 

aR 

Receiver  aperture  radius 

15  cm 

9R 

Receiver  f ield-of-view  1/2  angle 

2  mrad 

b 

Bistatic  receiver  offset 

0 

* 

aT 

Transmitter  aperture  radius 

3.1  cm 

9t 

Transmitter  divergence  1/2  angle 

1.2  mrad 

f 

Transmitter  focal  length 

-2583.33  cm 

• 

Eo 

Transmitter  pulse  energy 

-2 

10 

T 

Pulse  duration 

100  nsec 

• 

Medium 

/ 

a 

Backscatter-to-extinction  ratio 

.001 

e 

a 

Ambient  extinction 

0.2  •  10  ^  Cm  ^ 

t 

eo 

Peak  cloud  extinction 

3.4,  6.8,  17,  34  •  10~A  cm'1 

zo 

Cloud  center 

1  Km 

o 

Cloud  1/e  radius 

50  m 

.-„T~  v  “ .  .’•V. 
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3.  METHODS  OF  ESTIMATING  OBSCURANT  CLOUD  TRANSMISSIVITIES 
AND  CONCENTRATIONS  USING  LIDAR. 


3.1  Approach 

One  of  the  main  goals  of  the  present  contract  effort,  namely  the  de¬ 
velopment  of  a  mathematical  lidar  model  with  enough  flexibility  to  analyze 
a  broad  range  of  systems  and  obscuring  medium  conditions,  was  addressed 
in  the  first  part.  It  provides  a  framework  for  predicting  the  performance 
of  a  lidar  system  operating  against  three  dimensional  spatially  varying 
obscurants.  The  other  goal,  that  of  defining  a  viable  method  of  using  a 
lidar  to  infer  propagation  path  transmissivities  and  spatial  concentrations 
of  obscurants,  will  be  addressed  here. 

There  are  two  main  applications  for  the  data  on  obscuring  clouds  that 
can  be  derived  from  lidar.  One  use  is  to  estimate  the  transmissivity  be¬ 
tween  any  point  along  the  propagation  path  and  the  lidar  system.  This  is 
of  primary  interest  to  those  concerned  with  how  well  "observers"  either 
human  or  electro-optical  can  see  through  battlefield  obscurants.  The 
other  application  is  essentially  that  of  using  lidar  as  a  diagnostic  tool 
for  mapping  the  spatial  and  time  development  of  smoke  and  dust  generated 
by  explosive  rounds.  The  first  application  is  of  more  immediate  tactical 
interest  and  is  by  far  the  simpler  to  achieve. 

We  first  show  that  assuming  the  backscatter-to-extinction  ratio  a 
is  approximately  constant  leads  to  a  very  simple  means  of  estimating  the 
path  transmissivity  from  the  lidar  return  if  a  itself  can  be  estimated. 

The  sensitivity  of  the  proposed  approach  to  uncertainties  in  a  and  system 
noise  is  then  investigated,  and  a  generalization  of  the  approach  using  a 
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point  calibration  measurement  within  the  obscurant  cloud  is  developed 
that  would  permit  particle  concentrations  and  "CL"  values  to  be  determined 
from  the  lidar  as  well  as  transmissivities.  Finally,  limited  considera¬ 
tion  is  given  to  the  possibility  of  using  multiple  wavelength  systems  to 
enhance  the  reliability  and  range  of  usefulness  of  lidar  or  as  an  alter¬ 
native  to  the  point  calibration  method. 
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3.2  Derivation  of  the  Transmissivity  Algorithm 

Following  the  model  results  described  in  the  first  part,  an  approxi¬ 
mate  short  pulse  form  of  the  lidar  equation  for  a  Gaussian  aperture  and 
field-of-view  receiver  is 


Vt-2f) 


a  g(z)  r(z) 


aR  +  Qlz2+  <*T  +  9Tz)'2 


exp 


4  +  4  ^  +  (aT  +  6tz)2 


(27) 


where 


and 


A  s  t  *<-»  4  4  Eo  To 


T(z)  -  exp 


Z 

•/ 


dz‘  e(z’) 


0 


Eq.  27  results  from  Eq.  19  choosing 


(28) 


(29) 


E0  T0 

E(t)  -  -V2  5(t) 


(30) 


and  approximating  e(r)  “  e(e3z)  for  z  -  ct/2.  EQ  and  tq  represent  the 
pulse  energy  and  lidar  system  transmission.  This  equation  can  be  inte¬ 
grated  to  find  T(t)  is  we  assume  B(z)  u  a  e(z)  over  the  propagation  path; 
i.e.,  the  backscatter  is  uniformly  proportional  to  extinction.  Since 


e(z)  exp 


z 

.  z 

-2  I  dt*  e(z’) 

■  -1/2  h exp 

-2  /  dz*  e(z’) 

0 

0 

(31) 


integration  gives 
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3. 3  Sensitivity  Analysis  of  the  Transmissivity  Algorithm. 

Equation  32  allows  the  transmissivity  to  be  estimated  as  a  function 
of  propagation  range  through  the  obscuring  medium  as  an  indefinite  integral 
of  the  received  lidar  power.  To  make  the  idea  practical  several  things 
are  required: 

1)  It  is  necessary  to  have  an  accurate  knowledge  of  the  lidar  system 
parameters.  Note,  however,  the  absolute  value  of  the  power  collected  does 
not  need  to  be  known  because  of  the  normalization  A  which  includes  the 
transmitted  pulse  energy;  only  the  relative  power  received  to  that  trans¬ 
mitted  need  be  recorded. 

2)  Since  the  range  and  time  resolution  are  related  as  At  ■  2Az/c, 

100  ft.  range  resolution  requires  200  ns  resolution  in  the  recorded  lidar 
signal;  10  ft.  would  require  20  ns  resolution.  This  necessitates  a  short 
time  response  detector  and  high  bandwidth  recording  system. 

3)  Eq.  32  requires  an  estimate  of  the  backscatter-to-extinction 
ratio  a  be  made.  Since  this  is  often  poorly  known,  it  is  essential  to 
understand  the  effect  of  uncertainties  in  a  on  the  inferred  transmissivity. 

By  differentiating  Eq.  32  with  respect  to  a,  it  follows 


AT 
T  ' 


(33) 


Figure  4  plots  this  relation  for  various  values  of  T.  Clearly,  the  esti¬ 
mate  of  T  is  less  sensitive  to  errors  in  o  for  larger  T.  The  implication 
for  using  the  technique  on  dense  clouds  is  that  an  additional  local  measure¬ 
ment  of  a  may  be  necessary. 

In  order  to  explore  the  effects  of  various  uncertainties  on  the  esti¬ 
mates  of  T  using  Eq.  32,  a  computer  program  was  developed  to  compute  the 
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lidar  power  from  Eq.  27,  add  Gaussian  noise  to  simulate  real-world  con¬ 
ditions,  and  invert  the  signal  plus  noise  with  Eq.  32.  The  features  we 
wished  to  investigate  were  the  dependence  of  the  computed  transmissivity 
or.  noise,  cloud  density,  and  errors  in  a. 

Figure  5  plots  the  theoretical  lidar  signal  from  a  Gaussian  cloud 
superimposed  on  the  simulated  signal  with  Gaussian  noise  chosen  to  repre¬ 
sent  a  signal-to-noise  ratio  at  the  receiver  (z  m  0)  of  100  (top)  and  20 
(bottom).  The  corresponding  plots  of  theoretical  and  inferred  transmis¬ 
sivity  are  shown  to  the  right.  The  parameters  are  essentially  the  same 

-4  -1 

as  in  Table  1  with  Eq  ■  2.5  •  10  cm 

Examination  of  the  transmission  curves  shows  the  technique  is  capable 
of  inferring  the  correct  values  for  both  SNR's  up  to  the  cloud  center. 
Because  of  the  strong  signal  attenuation  from  returns  on  the  far  side  of 
the  cloud,  the  Inferred  transmission  estimates  are  unreliable  for  ranges 
greater  than  1  km.  Nevertheless,  detailed  inspection  of  the  numerical 
computer  output  shows  the  technique  can  produce  reasonably  correct  results 
for  signals  having  only  a  SNR  £  .1;  this  results  from  the  noise  averaging 
produced  by  the  integration. 

In  Figure  6  we  compare  the  inferred  transmission  for  two  peak 
Gaussian  cloud  densities  Eg  ■  10  ^  cm  ^  (top)  and  2.5  •  10  ^  cm  ^  (bottom). 
The  SNR  at  z  »  0  was  taken  to  be  100  in  both  cases.  Increasing  the  cloud 
density  by  a  factor  of  2,5  causes  drop  outs  in  the  inferred  transmission 
for  ranges  greater  than  the  cloud  center;  the  technique  works  best  for 
tenuous  clouds. 

Finally,  we  look  again  at  tne  problem  of  uncertainties  in  the  back- 
scatter-to-extinction  ratio  a.  Figure  7  shows  the  effect  of  producing 
simulated  lidar  signals  and  Inverting  them  with  10,  20,  and  50  percent 
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Sensitivity 


Sensitivity  of  lidar  inversion  to  cloud  density 


50%  error 


Sensitivity  of  lidar  inversion  to  uncertainty 
in  backseat ter-extinct ion  ratio. 
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errors  in  assumed  a.  For  these  runs  Eq  =  10  cm  and  SNR  *  100  at 
z  »  0.  There  is  comparatively  little  distortion  in  the  inferred  trans¬ 
mission  for  ranges  less  than  1  km  even  for  a  50%  error  in  a,  On  the  far 
side  of  the  cloud,  however,  a  10%  error  significantly  distorts  the  inferred 
transmission.  These  results  agree  with  earlier  observations. 

In  conclusion,  the  proposed  lidar  inversion  technique  should  be 
useful  for  relatively  weak  returns  (SNR  ~  .1)  if  the  backscatter-to- 
extinction  ratio  a  can  be  characterized  accurately.  Errors  in  a  will 
become  more  important  as  the  transmission  decreases. 
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3.4  Estimating  Obscurant  Concentrations  and  CL  Values. 

The  lidar  inversion  technique  described  previously  for  computing 
transmissivity  requires  an  accurate  estimate  of  the  backscatter-to- 
extinction  ratio  ot.  While  in  simple  cases  a  may  be  characterized  by 
a  priori  measurements  (e.g.  predominately  single  constituent  scattering 
from  a  known  material)  in  general  it  may  be  necessary  to  perform  a  local 
measurement  simultaneous  with  the  lidar  measurement.  If  this  is  possible, 
and  if  certain  other  parameters  can  be  measured  in  a  small  region  of  the 
propagation  path,  there  is  the  possibility  of  inferring  the  spatial  con¬ 
centration  throughout  the  region  of  non-zero  lidar  return.  The  purpose 
of  the  analysis  here  is  to  define  the  additional  local  measurements  needed 
to  calculate  the  spatial  concentration  and  CL  values  and  estimate  how  well 
the  teennique  would  be  expected  to  work  with  white  phosphorus  smoke  using 
either  a  1.06  pm  or  10.6  pm  lidar. 

From  Eq.  27,  the  extinction  coefficient  is  given  by 


v2f>  '4 +  4  *2  +  <*x +  v>2] 


A  a  T  (z) 


x  exp 


‘ _ £ _ 1 

+  0R  ^  +  (aT  +  9TZ)2J 


where  T(z)  is  estimated  from  Eq.  32.  To  compute  the  concentration  set 


Vz)  ■  p0 


Jdr  n(r,z)  »tot(r) 

0  00 

)  /dr  n(r,z)  (ff  r3) 


(36)  | 
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where  cr  fc(r)  is  the  extinction  cross  section  for  particles  of  radius  r 
and  n(r,z)  is  the  number  of  particles  of  radius  r  at  range  z  per  unit 
volume.  Here  is  the  density  of  the  smoke  or  dust  constituent,  assumed 
to  be  approximately  spherical  in  shape.  Set  N(z)  equal  to  the  total  number 
of  particles/volume: 


Then 


N(z) 


=  J dr  n(r , 


z) 


e(z)  -  N(z)  <{ o  y 


(37) 


(38) 


Cg(z)  -  pQ  N(z)<^Vp^> 


so 


(39) 


(40) 


where  and  are  the  average  extinction  cross  section  and  particle 

volume  for  the  density  function  n(r,z).  Finally, 


Equations  (40)  and  (41)  allow  CL  and  the  concentration  to  be  estimated 
using  Eqs.  (32)  and  (34). 

Applying  this  model  to  lidar  data  requires  estimates  to  be  made  in  a 
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localized  region  of  the  propagation  path  of  pQ,  a,  <^vp)>  and  <^tot^>  *  since 
the  material  of  the  assumed  single  scattering  component  is  presumably 
known,  Pq  is  known.  From  local  scattering  estimates,  the  bac. .scatter- 
extinction  ratio  a  can  be  estimated.  The  last  quantities  can  be  estimated 
from  local  measurements  of  n(r,zQ)  and  Table  2  summarizes  the 

steps  needed  to  estimate  C  (z)  and  CL  from  the  lidar  data. 

s 

For  the  analysis  here  we  used,  white  phosphorus  smoke  for  which  Pq  * 

1.3  g/cm3.  Using  a  mean  particle  radius  of  .57  pm  gives  /Vp^>  “ 

—3  ' 

.755  •  lO-3-3  cm3.  For  the  average  extinction  cross  section  we  used  10  cm 

and  2.94  •  10~8  cm3  at  10.6  pm  and  1.06  pm,  respectively.  The  smoke 

cloud  was  modeled  as  a  three-dimensional  Gaussian  density  centered  at 
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1  km  with  a  peak  number  density  of  3.4  *  10  particles/cm  and  a  one- 
sigma  of  50  m.  The  assumed  lidar  parameters  are  summarized  in  Table  3. 

Figure  8  plots  the  theoretical  and  inferred  CL  and  concentration 
values  for  the  two  different  wavelength  systems.  The  Inferred  values  were 
computed  from  the.  equations  above  using  simulated  lidar  returns  with  added 
Gaussian  noise  chosen  to  represent  a  NEP  of  1.372  •  10  watts.  (This 
works  out  to  a  SWR  of  103  at  the  transmitter-receiver.)  The  dropouts  in 
the  1.06  pm  plots  occur  at  points  where  the  inferred  transmission  becomes 
negative  due  to  the  simulated  noise.  It  is  clear  that  the  10.6  pm  system 
is  preferable  under  the  single  scatterer  conditions  assumed  here.  More 
work  is  needed  to  examine  the  advantages  of  using  both  wavelengths  for 
multiple  constituent  obscurants  with  broad  range  of  particle  sizes  since 
10.6  pm  return  may  be  too  weak  for  reliable  inversion  of  small  particle 
returns,  in  which  case  the  1.06  pm  lidar  would  be  preferable. 


TABLE  2.  LOGICAL  FLOW  OF  CALCULATION  OF  C  (z)  AND  CL  USING  LIDAR. 

g 


Known  scattering  constituent 
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SYMBOLS  USED  TO  INFER  CONCENTRATION  FROM  LIDAR 


Lidar  System 
Meaning 

Lidar  return  power 
Receiver  aperture  radius 
Receiver  field  of  view  h  angle 
Transmitter  beam  divergence  %  angle 
Transmitter  pulse  energy 
Transmitter/receiver  transmission 


Values 
Computed 
5.  cm 
.375  wr ad 
.5  mrad 
10-2  j 
.45 


Medium 


3.5  Multiple  Wavelength  Lidar 


As  stated  previously,  the  utility  of  the  proposed  lidar  inversion 
algorithm  for  estimating  transmissivities  and  spatial  concentrations  is 
greatly  enhanced  by  an  auxiliary  measurement  of  the  backscatter-to-extinc 
tion  ratio  a.  The  last  section  outlined  a  technique  incorporating  a 
local  measurement  of  scattering  quantities  including  a  to  normalize  the 
lidar  signatures.  For  those  cases  where  this  direct  measurement  is 
infeasible  some  other  method  must  be  employed  in  general  to  estimate  a. 

In  this  section  we  indicate  a  possible  approach  using  a  multiple  wave¬ 
length  lid*jr  system. 

If,  Instead  of  the  single  wavelength  lidar,  we  have  a  number  of 
simultaneous  wavelength  measurements,  the  backacatter  coefficient  for 
wavelength  X  at  range  z  becomes 


h(z)  " 


IX 


<2f)[a l  +  z2  +  (am  +  8^z) 


R  R 


T  T 


A  tJ(z) 


a 


exp 


2  2  2  ,  ,2 
aR  +  eR  1  +  ST  +  V> 


(42) 


where  is  computed  from  the  analogous  wavelength  dependent  version  of 
Eq.  32  and  an  approximate  estimate  of  ct. 

We  now  model  as 


sx<*> 


■  Jit  ntr. 


da 

z)  dir  (x’r’1t) 


(43) 


where  n(r,z)  is  the  number  of  particles  of  radius  r  at  range  z  per  unit 


volume  and  6.0 JAU  is  the  differential  backscattering  cross  section  for 
particles  of  radius  r  (assumed  spherical  here)  at  wavelength  X.  Equation 
43  constitutes  an  integral  equation  at  each  range  z  which  can  in  principle 
be  solved  for  n(r,z)  using  3^  from  Eq.  42. 

If  n(r,z)  can  be  estimated  from  Eq.  43  the  extinction  coefficient 
e(z)  can  be  computed  using  Eq.  35  and  estimates  of  atot(r)*  From  this 
the  transmissivity,  a,  concentration,  and  CL  values  can  be  computed  using 
the  equations  of  the  last  section. 

The  requirements  needed  to  solve  Eq.  43  for  n(r,z)  are: 

1)  The  kernel  of  the  integral  equation  do  /dft  must  be  specified 
using  either  Mie  theory  or  a  statistical  generalization  to  account  for 
different  particle  shapes  and  species. 

2)  Enough  wavelength  measurements  must  be  performed  with  high 
enough  slgnal-to -noise  ratio  to  provide  sufficient  degrees  of  freedom 

to  attempt  numerical  inversion  of  Eq.  43,  Both  of  these  are  very  strong 
requirements  and  more  study  is  necessary  to  determine  if  they  are  within 
the  grasp  of  current  theoretical  and  hardware  capability. 
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4.  SUMMARY  AND  CONCLUSIONS 


The  accomplishments  under  the  present  contract  effort  include 

1.  Development  of  a  three  dimensional  model  for  computing  the  lidar 
return  for  multiple  orders  of  scattering  for  arbitrary  spatially  varying 
media. 

2.  Computer  implementation  of  1.  for  the  case  of  first  order  scatter 
of  a  coherent  (Gaussian)  transmitter  beam  with  arbitrary  divergence,  a 
(Gaussian)  aperture  and  arbitrary  f ield-of-view  receiver,  and  a  diffuse 
reflecting  target  with  arbitrary  orientation  relation  to  the  incident 
radiation. 

3.  Derivation  and  sensitivity  analysis  of  a  means  of  computing  battle¬ 
field  transmission  by  integrating  the  collected  lidar  signal. 

4.  Analysis  of  a  means  for  estimating  spatial  concentrations  and  CL 
values  of  battlefield  obscurants  using  results  of  3.  and  a  local  measure¬ 
ment  of  scattering  parameters. 

The  motivation  for  constructing  the  lidar  model  was  to  provide  a 
predictive  tool  simple  enough  to  provide  physical  insight  into  the 
scattering  and  collection  process  yet  general  enough  to  treat  a  wide 
range  of  applications  and  model  extensions  without  reformulating  new 
equations  from  first  principles. 

The  Green’s  function  expansion  for  multiple  orders  of  scattering  of 
the  radiance  was  the  approach  chosen  since  it  incorporates  most  but  not 
all  of  the  possible  physical  effects  of  the  problem.  The  phenomena  it 
does  not  include  are  diffraction  and  polarization;  it  was  felt  that  these 
effects  are  relatively  unimportant  in  comparison  with  the  task  of  Including 
them  properly.  Emphasis  was  placed  on  the  first  order  return  since  for 
typical  smoke  or  dust  clouds  the  backseat ter-to-extinct ion  ratio  a  is 


39 


about  .001.  There  are  obscurants,  however,  which  could  conceivably  be 
of  tactical  interest  such  as  petroleum  oil  smoke  which  absorbs  relatively 
little  energy  in  the  ir*® and  would  require  more  than  the  first  order 
term  to  model  their  return.  In  this  case  the  general  formalism  presented 
for  higher  orders  would  be  of  more  than  academic  interest. 

The  lidar  model  was  applied  to  the  inverse  problem  of  inferring  in¬ 
formation  about  the  scattering  medium  from  the  lidar  return.  We  showed 
that  the  assumption  of  direct  proportionality  between  backscatter  and 
extinction  leads  to  a  simple  integration  of  the  first  order  lidar  equa¬ 
tion  to  give  the  path  transmission.  A  sensitivity  analysis  of  the  technique 
showed  it  to  be  relatively  insensitive  to  noise  but,  depending  on  the  cloud 
density,  quite  sensitive  to  the  assumed  backscatter-to-extinction  ratio. 
Under  the  assumption  that  a  localized  measurement  of  the  obscuring 
medium  parameters  could  be  made  simultaneous  to  the  global  lidar  measure¬ 
ment,  it  was  demonstrated  that  not  only  u  and  the  transmission  but  spatial 
concentrations  and  CL  values  could  be  measured. 

Finally,  we  examined  briefly  the  possibility  of  using  an  alternative 
to  the  local  measurement  involving  a  multiple  wavelength  lidar.  Clearly, 
it  is  preferable  to  perform  a  direct  measurement  of  the  scattering  para¬ 
meters  when  possible,  and  there  are  important  questions  remaining  to  be 
resolved  before  it  can  be  said  with  confidence  that  the  multiple  wave¬ 
length  approach  can  produce  comparable  results.  These  questions  involve 
modeling  the  differential  cross  sections  of  the  aerosols  constituting  the 
obscurant  and  examining  the  resulting  wavelength  sensitivity  to  see  how 

many  and  which  wavelength  are  required.  Furthermore,  multiple  scattering 

« 

components  introduce  substantial  additional  complications  which  must  be 
overcome  to  make  the  multiple  wavelength  approach  effective. 
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Appendix 


LIDAR  RETURN  FROM  FLAT  LAMBERTIAN  TARGET 


Following  the  approach  of  section  1.1  we  define  a  target  Green's 
function  in  the  small  angle  approximation  for  z  -  z*  «  0  as 


cT(£,n,t;£' ,n' ,t')  -  ~  xT(n  •  n^.)  BT(£.'  +  n'c )) 


x  T  6  (p  -  p'  -  (0  +  0')  c(t  -  t ' )/2)  A-l 


where  is  the  unit  vector  specifying  the  orientation  of  the  target  and 
6^  defines  the  target  reflectance.  The  transmission  T  is  given  by 

r  (t + 1 * )/2  i 

T  ■  exp  I -c/*  dr  e(p'  +  n'  c(t  -  t’))l  .  A-2 


For  a  Lambertian  target  the  scattering  function  x,p  is 


*T (n  ♦  n^)  - 


cos(n  •  n^) 


We  model  the  target  as  an  infinite  flat  surface.  The  condition  for 
to  be  non-zero  is  that  £_'  +  n'  c(t  -  t’)/2  lie  on  the  surface.  This 
is  expressed  mathematically  as 


[£.'  +  n*  c(t  -  t')/2  -  rT]  •  nT  -  0 


where  rT  ■  [0,  0,  z^]  is  tlie  vector  to  the  Intersection  of  the  plate  and 
the  z  axis.  Solving  for  £T  =  c(t  -  t') 
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ST  “  2^  -  p')  •  /<n'  •  n^) 


eT  -  eTo  <S(£t  -  c(t  -  t’))  . 


Proceeding  as  In  section  1,  we  define  a  transmitter  source  function 
using  Eq.  18  and  compute  the  target  radiance  at  the  receiver 


*  c  J "  dt  J"  d2£-'  •t’'^ 


where 


,  E(t  -  U c) 

2  XT(n  *  Et)  3t  T - 2 — 2 - 

0  it  Y 


P-  exp  [-(£  -  IV^/V2^]  A‘7 


£t  "  z[£-  “  It]  *  %  *  — t> 


y  -  1  -  5T/2f 


The  transmission  is  approximately 


T  *  exp 


r 

I  -  /*d£»  c(e3aT  -  £’>> 

"  Z_,  - 


For  the  target  oriented  perpendicular  to  the  transmitter  beam 
Jlj.  *  [0,  0,  -1],  £,j, ■  -2Zy/n^  a  2z,^.  Using  Gaussian  weighting  functions 
for  the  receiver  aperture  and  field-of-view 


t 


» 


I 


i 


ft 


ft 


ft 
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PT(t) 


2,  2 

"P  /aR 


A 


.-X 


c 

2 


—  T2 


E(t  -  2zt/c)  Tta^ 

au  +  (a  -  6  z  )2  +  e2z2 
Tt  ^  T  IT  R  T 


A-ll 


with 


T  ■  exp 


e(e3(2zT 


5*)) 


A-12 


For  non-zero  target  orientations,  the  integration  of  Eq.  A-ll  is  much 
more  difficult.  The  case  of  a  square  pulse  can  be  integrated  in  closed 
form  in  terms  of  error  integrals  but  more  general  cases  require  machine 
integration.  Figure  2  plots  the  result  for  a  20  ns  square  pulse  with 
various  target  orientations. 
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ft 
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